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Abstract 

We construct approximate solutions to the time-dependent Schrodinger equation 

for small values of fi. If V satisfies appropriate analyticity and growth hypotheses 
and \t\ < T, these solutions agree with exact solutions up to errors whose norms are 
bounded by 

C exp { - j/H } , 

for some C and 7 > 0. Under more restrictive hypotheses, we prove that for sufficiently 
small T', \t\ < T" \ \og(h)\ implies the norms of the errors are bounded by 

C exp{- 7 7r}, 

for some C, 7' > 0, and a > 0. 
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1 Introduction 



In this paper, we construct exponentially accurate semiclassical approximations ip(x, t, h) to 
certain normalized exact solutions *&(x, t, K) of the (i-dimensional time-dependent Schrod- 
inger equation 

h 2 

lh°— = -yA$ + V*. (1.1) 

More precisely, our main result is that for \t\ <T and small values of h, these approximations 
satisfy error estimates of the form 

|| ^{x, t, h) - *{x, t, h) \\ L2{Rd) < Cexp{- 7 /^}, (1.2) 

where 7 > 0. 

Our construction of ip(x, t, h) is technically complicated, but quite explicit. It uses a 
particular collection of semiclassical wave packets {(pj(A, B, h, a, r], ■)} that are defined 
in |J, J7J, and the next section. Here A and B are d x d complex matrices that satisfy 
certain conditions. The quantities a and 77 are elements of IR d . For fixed A, B, Ti, a, and rj, 
{ tfj(A, B, h, a, rj, • ) } is an orthonormal basis of L 2 (IR d ) as j ranges over all <i-dimensional 
multi-indices. 

The function fj(A, B, %, a, rj, • ) is concentrated near position a, and its Fourier trans- 
form is concentrated near momentum 7]. Its position and momentum uncertainties are pro- 
portional to V%. The position uncertainty is determined by the matrix \A\ = \/AA*, and 
the momentum uncertainty is determined by \B\ = yBB* . 

We construct ip(x, t, h) by applying the idea of "optimal truncation" of an asymptotic 
expansion. For initial conditions of the form 

0, ft) = £ Cj ^(A(0),B(0),h,a(0), V (0),x), (1.3) 

\j\<J 

with J2\j\<j \ c j\ 2 = 1> there exist (f5|, ||) approximate solutions 

^(x, t, ft) = e iS W h J2 c i(l, ft) <Pj(A(t), Bit), ft, a{t), r,(t), x), (1.4) 

mm 

that satisfy 

sup \\^{x, t, h) - t, h)\\ L 2 {Ri) < C{1) h 1 ' 2 (1.5) 
te[-T,T] 

for some constant C(l). Here J (I) = J + 3/ — 3, and A(t), Bit), ait), rj{t), and S(t) are 
solutions to the classical equations of motion 

= 

17(f) = -W(a(f)), 

Ait) = iB{t), (1.6) 

Bit) = iV {2 \ait))Ait), 

S(t) = ^-Via(t)), 
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where denotes the Hessian matrix for V, and the initial conditions A(0), B(0), a(0), 
77(0), and 5(0) = satisfy 

A t (0)B(0) - B t (0)A(0) = 0, (1.7) 
A*(0)B(0) + B*(0) A(0) = 21. (1.8) 

The Cj(l, t, h) satisfy a linear system of ordinary differential equations that we describe in 
the next section. 

We carefully estimate the /-dependence of C(l) in (|1.5|) . Then, for each h, we choose 1(h) 
to minimize the error C(i(h))h m/2 over all choices of I. It turns out that 1(h) behaves like 
a constant times 1/h. We define ip(x, t, h) = ipi^(x, t, h) and prove that (|1.2|) is satisfied. 

For t in a fixed compact interval, the precise statement of our results is the following: 



Theorem 1.1 Suppose V is a real-valued function on IR d that is bounded below and has an 
analytic continuation to the set 

D = {z E(U d : \Im zA < S, j = 1, 2, • • • , d}. (1.9) 

Suppose further that there exist M > and t > 0, such that 

\V(z)\ < M exp(r\z\ 2 ), for all z G D, (1.10) 

where \z\ 2 = J2j=i \ z j\ 2 ■ Suppose initial conditions A(0), B(0), a(0), r)(0), 5(0) = ; and 
c j f or \j\ — J are specified that satisfy (\1-Tj) and J2\j\<j \ c j\ 2 — 1- Then for any T > 0, 
there exist C and 7 > 0, such that the difference between the semiclassical approximation 
ip(x, t, h) and the exact solution ty(x, t, h) to the Schrodinger equation ( |i.i[ ) with initial 
condition satisfies 

|| ip(x, t, h) - ty(x, t, h) || L 2 (Rd) < C exp{-7//i}, 

whenever \t\ < T. 

Remarks 1. Theorem |1.1| can be generalized to allow time-dependent potentials. For 
example, suppose a potential V(x, t) depends smoothly on t, is bounded below, and satisfies 

±\D»V(x,t))\ < Mex P^I 2 ) 



5\m\ 



for \t\ < T for all multi-indices m. Suppose further that a classical solution (a(t), r)(t)) to 
Newton's equations with potential V(x, t) is bounded for \t\ < T. Then the conclusion to 



Theorem 1.1 holds. 



2. Our results can also be extended to obtain weaker error estimates of the form 
C exp { — r j/h tJ } for some a G (0, 1), when the potential belongs to a Gevrey class. 



3. Theorem |1.1| is optimal in the sense that the conclusion fails if the hypotheses are relaxed 
slightly. For example, consider the one-dimensional potential V(x) = exp(— l/x u ), for x > 0, 
and V(x) = for x < 0, where u > 0. It is shown in M that this potential belongs to the 
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Gevrey class of order 1 + 1/u. For initial conditions a(0) = 0, 77(0) = 0, ^4(0) = -B(O) = 1, 
S(0) = 0, and Cj(0) = 5jfl, our approximation yields (fo(l + it, 1, h, 0, 0), for all times. This 
function is very simple, and we can write the error term explicitly. By steepest descent 
analysis, we can show that there exist 5 > and Si > £2 > 0, such that t G (0, h s ) implies 

exp(-S 1 /^/( 1+ ")) < || e- itH ^' h <p (l, 1, h, 0, 0, ■) - tp (l + it, 1, % 0, 0, •) || 

< exp(-£ 2 /ft u/(1+u) ). 

Note also that if we choose a(0) = —a < and 77(0) = rj > 0, it is easy to check that the 
error term is 0(exp(— r y/K)), for each t < a/r). 

4. For all practical purposes, we can replace the c/s by the corresponding Dyson expan- 
sion up to order 1(h) without spoiling our exponential estimate. The normalization of the 
approximation, however, will be lost. 



Theorem |1.1| can also be generalized to allow time intervals that grow like | \og(h)\ as h 
tends to zero. However, we obtain a somewhat weaker conclusion. Our precise results are 
summarized by the following theorem. 



Theorem 1.2 Suppose V is bounded below and analytic in 

D = {z e(U d : \Im Zj\ < 5, j = 1, 2, d}. (1.11) 
Suppose further that there exist M > and t > 0, such that 

\V(z) I < M exp(r|z|), for all z E D. (1.12) 

Suppose initial conditions A(0), B(0), a(0), 77(0) , S(0) = 0, and Cj for \j\ < J are specified 
that satisfy (\1.7\) and J2\j\<j \ c j\ 2 — 1? and further assume there exist N > andX > 0, such 
that 

\\A(t)\\ < Nexp(X\t\). (1.13) 

Then for sufficiently small X" > ; there exist C , 7' > ; and o > 0, such that the difference 
between the semiclassical approximation ip(x, t, h) and the exact solution \l/(x, t, h) to the 
Schrddinger equation ( |i.i[ ) with initial condition satisfies 

|| 1p(x, t, K) — ^(x, t, h) \\L 2 (R d ) 

whenever \t\ < T'\log(H)\. 

Remark Standard existence and uniqueness theorems for systems of ODE's show that 
condition ( |1.13|) is satisfied if the norm of the Hessian V^(a(t)) is uniformly bounded. That 



is the case if V is the sum of a quadratic polynomial plus a bounded analytic function. It is 
also the case if E denotes the energy of the considered trajectory and ||V^ 2 )(x)|| is bounded 
on the connected component of the classically allowed region De = {x G IR d : V(x) < E} 
that contains a(t). This is satisfied for all confining potentials. 
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The propagation of coherent states is also considered by Combescure and Robert in [[|, 
using an approximation given by a linear combination of squeezed states. (The squeezed 
states coincide with our semiclassical wave packets, although the notation is quite different.) 
Their emphasis is on the long time behavior of this approximation. The bound on the error 
term is of the form Ci(t)h 1 ^ 2 , with expicit control of the time-dependence of Ci(t) in terms 
of classical quantities. The I behavior is however not investigated. 

Results of a flavor similar to ours can be found in the work of Yajima fllTf . They are 
obtained by means of the pseudo-differential techniques developed in the analytic context 
by Sjostrand in flQ| . These results concern the propagation of wave packets of the form 
ip(x) = e tS ( x )/ h f(x), where S is analytic and / belongs to the set of compactly supported 
Gevrey functions of order s > 1. Assuming the potential V is analytic, Yajima constructs 
approximations to the actual evolution of such wave packets that are valid up to an error, 
whose L 2 (IR d ) norm is of order e~'y/ hl/i2a 1) ; with 7 > (see Theorems 1.2, 1.2 and Lemma 
2.5 in |11[| ). However, it should be possible to make use of the theory jnj to recover our 
results. 

Similar issues have been dealt with by Bambusi, Gram and Paul in 0. They focus on 
the validity for large times of the semiclassical approximation of the Heisenberg evolution of 
a smooth observable, under analyticity assumptions on the hamiltonian. They prove that 
the semiclassical approximation remains useful for times up to order | ln(ft)|, the Ehrensfest 
time scale. However, the Hamiltonians they can accomodate consist more or less of analytic 
perturbations of the harmonic oscillator that decay as x and p tend to infinity. 



The paper is organized as follows: In the Section 2, we prove Theorem |1.1| under the 
assumption that two types of error terms satisfy certain bounds. We prove the two required 



bounds in Sections 3 and 4. In Section 5 we describe the proof of Theorem 1.2 



2 Proof of Theorem 1.1 

We begin this section by presenting the definition of the semiclassical wave packets 

<Pj(A, B, h, a, f], x) that is given in 0. A more explicit, but more complicated definition is 

given in ||. Since Jf| provides a detailed discussion of these wave packets, we do not prove 

all their properties here. 

We adopt the standard multi-index notation. A multi-index j = j 2 , . . . , jd) is a 

d-tuple of non-negative integers. We define \j\ = ELi h, 3- = O'lOC^O ■ ■ ■ O'dO, 

■ „• O 3 

x 3 = X1X2 ■ ■ ■ x J d d i and D 3 = — — — — — — — — . 

d (dxi)K(dx 2 ) 32 ■ ■ ■ (dx d ) 3 d 

Throughout the paper we assume a G IR d , rj e IR d and h > 0. We also assume that A 

and B are d x d complex invertible matrices that satisfy 

A l B - B l A = 0, (2.1) 
A* B + B* A = 21. (2.2) 

These conditions guarantee that both the real and imaginary parts of BA^ 1 are symmetric. 
Furthermore, KeBA -1 is strictly positive definite, and (KeBA^ 1 ) 1 = A A*. We note that 
conditions fl2.1|) and (|2.2| ) are preserved under the dynamics generated by ( |1.6|) . 
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Our definition of <fj(A, B, %, a, 77, x) is based on the following raising operators that are 
defined for m = 1, 2, . . . , d. 



r d 



Am(A,B,%,a,r))* 



2% 







^ ] B n m [p^n &n) ^ ] A nm ( %h V'n) 

Ln=l n=l c * r ™ 



Definition For the multi-index j = 0, we define the normalized complex Gaussian wave 
packet (modulo the sign of a square root) by 

<p (A, B, %, a, 77, x) = 7T- d/4 h- d/4 (det(A))" 1/2 

x exp{-((x-a), B A' 1 (x - a) }/(2h) + % (77, (x - a) )/%}. (2.3) 

Then, for any non-zero multi- index j, we define 

^(A,B,h,a, V ,-) = (A(A5,M,l)f (^A^.^M)? ••• 

x (A^A-BiM,?/ 1 )*)^ Po(4,.B, ft, a, 77, ■)■ 

Remarks 1. For A = B = I,h = l, and a = 77 = 0, the <fij{A, B, %, a, 77, • ) are just the 
standard Harmonic oscillator eigenstates with energies \j\ + d/2. 

2. For each A, B, %, a, and 77, the set { <Pj(A, B, %, a, 77, • ) } is an orthonormal basis for 
L 2 (iR d ). 

3. The raising operators can also be given by another formula that was omitted from [0] in 
the multi-dimensional case. If we set 

g(A, B, %, a, x) = exp { - ( (x - a), (EA -1 )* (x - a))/(2h) - i (77, (x-a))/%\, 

then we have 

[% 1 d — d 

( A m (A, B, %, a, 77)* ip)(x) = - \ - r A nm ^— ( g(A, B, h, a, x) ip(x) ) . 

V 2 g(A, B, %, a, x) ~ dx n 

4. In ||, the state <fj{A, B, %, a, 77, x) is defined as a normalization factor times 

Hj(A; h~ 1/2 (A)" 1 (x - a)) (p (A, B, %, a, 77, x). 

Here Tlj(A; y) is a |j| th order polynomial in y that depends on A only through Ua, where 
A = \A\ Ua is the polar decomposition of A. 

5. By the argument on page 370 of [|] or by scaling out the \A\ and % dependence and 
using Remark 3 above, one can show that Tij(A; y) e~ y I 2 is an (unnormalized) eigenstate of 
the usual Harmonic oscillator with energy \j\ + d/2. 

6. When the dimension d is 1, the position and momentum uncertainties of the 



<Pj(A, B, %, a, 77, • ) are J (j + l/2)h \A\ and J (j + 1/2)% \B\, respectively. In higher dimen- 



sions, they are bounded by J (|j| + d/2)% \\A\\ and J + d/2)% \\B\\, respectively. 



7. When we approximately solve the Schrodinger equation, the choice of the sign of the 
square root in the definition of ifo(A, B, h, a, rj, • ) is determined by continuity in t after an 
arbitrary initial choice. 

The proof of the theorem depends on the following abstract lemma. 



Lemma 2.1 Suppose H(h) is a family of self- adjoint operators for h > 0. Suppose ip(t, ti) 
belongs to the domain of H{Ti), is continuously differentiable in t, and approximately solves 
the Schrodinger equation in the sense that 



. . dip 

ih— {t, h) 



H(h)ip(t, h) + f(t, h), 

where h) satisfies 

U(t, h)\\ < n(t, %). 

Then, for t > 0, 

\\ e - itH{n)/h ip{o, h) - ip(t,h)\\ < h- 1 

The analogous statement holds for t < 0. 



/i(s, h) ds. 



(2.4) 



(2.5) 



(2.6) 



Proof: Assume t > 0; the proof for t < is similar. By the unitarity of the propagator 
e -itH(h)/n an( j fundamental theorem of calculus, the quantity on the left-hand side of 
TM can be estimated as follows: 



-itH(h)/h 



^(o, h) - h) 



V>(0, h) - e itH(n)/n ij)(t, h) 
* d 



, fV(o, h) - e lsH{h)/h if>(s, h)) ds 
o as v 7 

-i h- 1 e lsH(h)/h H(K) if>(s, h) - e 

f ih- 1 e isH{h)/n £(s, h) ds || 
Jo 



< h 1 / /j,(s,h)ds. 
Jo 



This proves the lemma. I 

Because V is smooth and bounded below, there exist global solutions to the first two 



equations of the system (|1.6|) for any initial condition. It then follows immediately that the 
remaining three equations of the system ( p..6| ) have global solutions. Furthermore, it is not 
difficult (0, 0) to prove that ( [2.1|) and (|2.2j ) are preserved by the flow. 
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As mentioned in the introduction, it is proved in |J and || that initial conditions of the 
form (|1.3|) give rise to approximate solutions of the form 

\i\<J(i) 

with errors whose norms are of order h l/2 . Here J (I) = J + 3/-3, and A(t), B(t), a(t), rj(t), 
and S(t) satisfy (|1.6|) . The coefficients Cj(l, t, K) satisfy the linear system 

ihc k (l,t,h)= K kj (i,t,h) Cj (l,t,h), \k\ = Q, l, J(l), (2.7) 

\j\<J(l) 

with initial conditions Cj(l, 0, h) = Cj, for \j\ < J in accordance with ( p.. 3D and Cj(l, 0, %) = 
for \j\ > J. 

To specify the ( J + 3/ — 2) x ( J + 3/ — 2) matrix K(l, t, h) that appears in (|2.7| ), we first 
decompose the potential as 

V(x) = W a (x) + Z a (x) = W a (x) + (V(x) - W a (x)), (2.8) 

where W a (x) denotes the second order Taylor approximation (with the obvious abuse of 
notation) 

W a (x) = V(a) + V {1) (a) (x — a) + V^(a)(x-a) 2 /2, (2.9) 

Next (reverting to multi-index notation), we approximate Z a (x) by its Taylor approximation 
of order / + 1, 

ZM(*)= £ ^^-(x-ar, (2.10) 

3<\m\<l+l m ' 

and define the infinite matrix 

K kj (l, t, h) = 

( <p k (A(t), B(t), h, a(t), n{t), x), Zf {t) {x) Vj {A(t), B{t), h, a{t),r](t),x) ). (2.11) 

Then, we obtain the matrix K(l, t, Ti) from K(l, t, Ti) by restricting the indices to \ j\ < J (I) 
and \k\ < J(l). 

The general strategy ||, ||, to show that ipi(x, t, h) is an approximation to the actual 
solution t, h) of and (|1.3| ) up to order ft^ 2 is as follows: From [fjj, we know that 
for all multi- indices j, 



■ft 9 



* ff W/* V » i (A(t),S(t),ft l o(t),»7(t),x) 



-y A + W^a;)) [e <s W/ R ^(^(*),5(t),^a(t),»7(t),x)] . (2.12) 
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Thus, the (pj take into account the kinetic energy and W a (t)(x) parts of the Hamiltonian. 
Next, we expand the exact solution as 

*(x, t,h) = J2 bj(h, t)e lS ^ h ^(A(t), Bit), h, a(t), V (t), x), 
j 

where the bj(h, t) satisfy an infinite linear system of ordinary differential equations whose 
matrix is obtained from the Z a ^ (x) term in the Hamiltonian. In that system, we make a first 

approximation by replacing the function Z a ^(x) by its Taylor approximation Z^j+Jx). This 

yields an infinite linear system whose matrix is K(l, t, h). Its entries are time dependent 
polynomials in h 1 ^ 2 of order I — 1. We make a second approximation by truncating the 
infinite system to obtain (|2.7|) that is satisfied by the Cj(l, t, h). The result (|1.5|) is proved 
by using Lemma ^TT] to show that the errors generated by the Taylor approximation and the 
truncation approximation are of order h 1 ^ 2 . 

As described in the introduction, we construct the exponentially accurate approximate 
solution i/j(x, t, h) = ?pi(h)(x, t, h), by keeping track of the /-dependence of C(l) in ( |1.5| ) 
and then choosing 1(h) in order to minimize the error. 



In the remainder of this section, we prove Theorem 1.1 under the assumption that the 



following three technical lemmas are true. We prove the first lemma in Section 3. It estimates 
errors that arise from our replacement of Z a ^)(x) by Z^Ljx). The second and third lemmas 
are proved in Section 4. They bound certain matrix elements and combinatorial quantities 
that arise from the truncation approximation discussed above. 



Lemma 2.2 Suppose V satisfies the hypotheses of Theorem \l.\ \m\ = 1+2, J{1) = J+31—3, 
and 



4>i(x, t, h) 



JS(t)/h 



]T Cj (l, t, ft)^(A(t),5(t),M0O,^), 



X 



\j\<J(l) 



with Ey^jp) fail, t, h)\ 2 = 1. 

Let ((x, a(t)) = (a(t) + 9 x>a {t){x — a(t))) e IR d , with 9 X)0i ( t \ G (0, 1). There exist constants go 
and gi, that depend on d and J only, such that for sufficiently small h, 



D m V(C(x,a(t))) 



(x-a(t)) m ^(;t, h) 



< g M exp(4r(5 2 d + a(t) 2 )) ( 9l yJ%(l + 2) \\A(t)\\/5 



1+2 



(2.13) 



Lemma 2.3 We define the infinite matrix 
Xf k it,h) = 

( V] {A{t), B(t), H, a(t), V (t), •), (x - a(t)) m ^ k {A{t), B(t), h, a(t), V (t), ■) ), (2.14) 
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and then define the finite matrix X m (l, t, h) from X m (t, h) by restricting its indices to 



\j\ < J(l) and \k\ < J{1). Then, Xf h (l, s, h) = and X™ k (s, h) = if 
for each N > 0, there exists D(N) < oo, such that 



sup 

\k\ < J(i) 
J(i) + ± < b'l < J(i) + i + i 



(x m ° (t , h)x m > (i, t x , h)x m * (i, t 2 ,h)--- x m « (i, t q , h) 



> \m\ and, 



|mo| + |mi|+|m2|---+|m g | 



< \D(N)y/hJ(l) sup \\A(t)\\ 
\ te{to,ti,t 2 ,--tq} 



(2.15) 



for any collection tuq, mi, . . . , m q of multi-indices that satisfy \rrij\/J(l) < N . 



Lemma 2.4 We define 

F p (n, q) = Y, 1 ( 2 - 16 ) 

1 < \m\\, |rrea|, — , \m g \ < p 
\mi | + [ma | + ■ ■ ■ + \m g \ = n 

to be the number of distinct sets {mi, m>2, m q }, where each rrij is a d-dimensional 
multi-index with 1 < \rrij\ < p and \mx\ + |m 2 | + • • • + \m g \ = n. We note that F p (n, q) is 



zero un 



less q < n < qp. Suppose that a function L q (l) satisfies 



L ^)<W^ E (C 2 hir' 2 F l+1 (n, q), (2.17) 

11 <?• n=3l-2 

where C\ and Gi are constants. Let [a] denote the greatest integer less than or equal to a. 
Then there exists g* > 0, such that for any g e (0, g*), there exist positive constants C%, 71, 
and H* , that depend only on g, C\, and C%, such that 

0<h<h*, 1(H) = {g/hj, and 2 < q < 1(h) + 2 (2.18) 

imply 

L q (l(h)) < C 3 exp{- 7l /^}. (2.19) 

Proof of Theorem |1.1| : We define ipi(x, t, h) by (|1.4j), where the Cj(l, t, h) are determined 
by the system ( |2.7|) and initial conditions described above. 
To apply Lemma [2.1| we define 

h(x, t,h) = ih —^(x, t,h) - I - — A + V(x) j ^(x, t, h). (2.20) 

By using ( |2.12j ), we see that this can be decomposed as a sum of two terms, 

^\x, t, h) = (Z%(x) - Z a(t) (x)) fr(x, t, h) (2.21) 
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and 



(2.22) 



where -P{|j|>«+2} is the orthogonal projection onto the span of the set 
B{t), h, a(t), V {t),-): \j\ >l + 2}. 
By the standard Taylor series error formula, 

*<«<*)-*&■<»)= £ °" v( «T- a(t))) (»-.(«))■ 



m|=i+2 



for some ((x, a) = (a + 6 X)a (x — a)), with G (0, 1). Thus, by the crude estimate 

E i < E i < E E i 

m|=n |m|<ra J=l m,j<n 



[n 



+ l) d = e dln{n+1) < (e d y 



(2.23) 



and Lemma |2.2| , we obtain 

lk, (1) (-,*, *) 



< E 

|m|=Z+2 



D m K(C(x,a(t))) 



ml 



(x-a(t)) m 



< g M exp(Ar(6 2 d + a(t) 2 )) (g iy Jh(l + 2) \\A(t)\\/S 



1+2 V"' 



Thus, there exist constants C4 and C5, such that |i| < T implies 

|e, (1) o,*, &)| < c 4 (c 5V / m/ + 2)) z+2 . 

The quantity ^ 2 \-, t, %) satisfies 

{ip^A{f),B{t),n,a{t),ri{f),x), tf\x, t,h)) = 0, if 



(2.24) 



< \j\< J (I) or 
\j\>J{l)+l + l, 



and 



(<Pi(A(t),B(t),h,a(t),7i(t),x), tf\x, t, h)) = 
£ ( <Pj{A{t), Bit), h, a(t), r)(t), x), Z a(t) {x) <f k (A(t), B(t), h, a(t), r)(t), x) ) c k (t) 

\k\<J(l) 

= (K(t) c(t))j, if J(l) <\j\< J(l) + 1 + 1, (2.25) 

where we denote the Cj(l, t, %) collectively by the vector c(Z, t, Ti). We easily verify these 
facts by using (|2.7|) , (|2.11| ), (|2.12| ) and Lemma |2~3| . To estimate the norm of ^ 2 \-, t, h), we 
use the Dyson expansion with remainder to decompose 



c(Z, t,h) = J2 cq ( l , t, h) + r(Z, t, h), 

q=0 



(2.26) 
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where (dropping some arguments) 



C*(f) 



dsi / ds 2 



lg-1 



ds q {ih)-' 1 K{s 1 )K{s 2 )-- -K{s q ) c(0), 



(2.27) 



and 



r(t) = r^ 2 --- f Sl ds l+1 (m)-^K( Sl )K(s 2 )---K(s l+1 )c(s l+1 ). (2.28) 

Using (PUD , (127L0D , and (gTp , we see that each c q (l, t, %) is of order h q/2 and that r(l, t, %) 
is of order /i (z+1)/2 . 

To estimate the norm of ft) we study the j th component of K(t)c(l, t, h), with 

\j\ > J (I). Because of ( p. 26 ), this coefficient is a sum of two types of terms: those that arise 
from c q (l, t, h) and those that arise from r(Z, t, Ti). Using fl2.10| ), we expand K(l, t, h) in 
(P3) and (|^D and K(l, t, K) in terms of X m (l, t, Ti) and X m (t, h), to obtain 



_ p+i)(«+i) 

iif(t)c«(o = (iny q E 

n=3(g+l) 



E 



mo, mi , m2, ■ • • , m 9 
m | + |mi| + • ■ ■ + |m s | = n 
3 < |mj[ < J + l 



81 



dsi / d«2 ■ ■ 



and 



rfs £> m °y(a(t))£) roi y(a(g 1 )) ■ ■ ■ 73 m «1/(a(s ? )) 
x X mo (t)X mi ( Sl )X m2 (s 2 ) ■ ■ -X^isg) c(0), 



_ 0+1)0+2) 

K(t)v(t) = (ih)-<w E E 

n=3(i+2) mo,roi,m2r"i m I+i 

|m | + |mi| H + |m i+1 | = n 

3 < Kj| < £ + 1 



(2.29) 



«1 



dSi / C?S2 " " 



D mo V{a{t))D mi V{a{s 1 )) ■ ■ ■ J D m '+ 1 1/(a(s m )) 
dsi + i j j j j 

x X mo (t)X m ^ Sl )X m2 (s 2 ) ■ ■ -X m ^(s l+1 ) c(s,+i). 
The values of rrij that occur in both Q2.29Q and (|2.30|) satisfy 

\mj\/J{l) < (/ + l)/(3/-3) < 1, 
as long as I > 3. So, we can apply Lemma |2l| with N — 1. 



(2.30) 
(2.31) 



IjH*I 



> |m|. Since Cfc(0) is non-zero 



Recall that X^ k (l, s, h) = and Xf fc (s, a) = if 

only for \k\ < J, and we need only consider the j th coefficient of K(t)c(t) for \j\ > J(l), the 
only relevant values of n in (|2.29p must satisfy 

n>3Z-2. (2.32) 

This condition is also satisfied for all values of n in ( |2.30| ), since the sum begins with n = 3Z+6. 



12 



To use the analyticity assumptions to get estimates on the derivatives of V, we define 
C s (x) = {ze(U d : Zj = Xj + be ie \ Q j G [0, 2vr), j = 1, 2, ■ • • , d}. 
If z G C$(a(t)), then, for all j — 1, 2, • • • , d, 

\zj\ < 5 + \a,j(t)\. 



Hence, writing — - D m V(a{t)) as a <i-dimensional Cauchy integral, we get the bound 
ml 

±\D"V(a{t))\ < ^^ V ^ - «(*) 



ml 



s\™\ 



(2.33) 



where 



v(t) < M exp(4r(5 2 d + a 2 (t))). 
Furthermore, since ||A(£)|| depends continuously on t, there exists w(T), such that 

sup \\A(t)\\ < w(T). 

te[-T,T] 

The norm of the vector c(i) is 1 since K(l, t, K) is self-adjoint. Thus, the non-zero entries 
of c(0) are each bounded by 1, and by a crude estimate, there are at most (J + l) d of them. 
Similarly, c(i) has at most (J(l) + l) d non-zero entries, each of which is bounded by 1. Thus, 
for I > 3, (FTTID , (PT0|) , ( CT ), (CT), (|2730D , ( ggg) , (^3TD , and Lemma |p| imply the 
following two estimates when j satisfies J(l) + 1 < \j\ < J(l) +1 + 1: 



< v(t) 



(j T v(s)ds) 9 C+i)C«+U 



£ (D(l) w{T)/8) n h n ' 2 J{l) n ' 2 F l+1 (n, q + 1) ( J + If 



n=3l~2 



v(t)h(q + 1)(J + l) d (/o T ^(^)^) <?+1 



Jo v(s)ds 



. £ P(i)w(r)VM0/*) B ^Ui(n,(7 + i) 

19 + L J- n=3l-2 

(2.34) 



and 



(K(t)r(t)) 



< 



v(t)h(l + 2)(J(l) + l) d (lo v(s)ds) l+2 



Jo v(s)ds 



11 \l + 4)<- n =3l-2 

(2.35) 



where F p (n,q) is defined by ( |2.16| ). 
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By Lemma ( gig ), and ( ggg ), both c«(*)) 3 



and 



(*(*M*))j 



C3 exp { —71/^} for an appropriate choice of = [p/ft] and sufficient 



are bounded by 
y small H. Each 



of the I + 2 terms in (|2.26|) contributes a term of this type to Q (x, t, ft), so 



Ti 



-1 



^(x, t, n) < C 3 h- 1 (l(h) + 2) exp {-71A} 
< C 3 exp { -72/^} , 



for any 72 < 71 when fi is sufficiently small. 
We shrink g if necessary to make 

C\g<\ 

and set I = 1(h) in (|2.24|) . This yields a similar estimate 



(2.36) 



jP(;t, %)\\ < C 6 exp{-7 3 /^}, 
for some 73 > 0. 

We combine these two estimates and apply Lemma |2.1| to obtain (|1.2|) with 
7 = min{72, 73}. This proves the theorem. I 



3 Proof of Lemma 2.2 



For simplicity, we drop the t dependence in the notation throughout this section. 

To prove Lemma ^2| , we use Holder's inequality to see that J2\j\<Jm \cj(l, Ti)\ 2 = 1 
implies 





( Y /2 




E \cj(!,*)\ < 


E 1^, 


(?■) 


u\<j(i) 


\ \j\<j(i) J 


V \i\<Ai) ) 



1/2 



< (J{l) + l 

Thus, it is sufficient to prove 
D m V{((x,a)) 



d/2 



nr. 



(x — a) m <fj(A, B, h, a, 77, x) 



< g 3 M exp(4r(5 2 d + a 2 )) ( g 4 yjh{l + 2) \\A\\/S 



1+2 



(3.1) 



for some g 3 and g±, whenever \j\ < J (I). 

We mimic the proof of |2.33| to obtain a bound on D m V(((x,a))/m\. If z G Cs(((x,a)), 
then, for all j = 1, 2, • • • , d, 

\zj\ < 8 + \Q(x,a)\ < 6 + \a,j\ + \xj — a,j\. 

Using this and applying (b+c) 2 < 2(b 2 + c 2 ) several times, we see that z G Cs(((x, a)) implies 

| V{z) | < M exp(2r(a; - a) 2 ) exp(4r(5 2 rf + a 2 )). 
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Hence, writing — - D m V(((x, a)) as a d-dimensional Cauchy integral, we easily obtain the 
ml 

bound 

jL \D™V(((x, a))\<M 6Xp(4r( ^ + Q2)) exp(2r(x - a) 2 ). (3.2) 



Thus, estimate (|3.1| ) follows from the corresponding estimate on the integral 

(a; - a) 2m exp(4r(x - a) 2 ) {cp^A, B, h, a, rj, x)\ 2 dx. (3.3) 



IIR d 

Performing the change of variables x i— > ?/ = l^l-^x-a)/^ 1 / 2 , and using the explicit formula 
for <fj, we see that 

1 = abij!^ y^(l A ^) 2mex P(^ 2 + 4 ^(l A ^) 2 ) |Wi(^;y)l a dy, (3-4) 

where T~Ci(A; y) is the polynomial described in Remarks 4 and 5 that immediately follow the 
definition of <fj{A, B, H, a, rj, x) in Section 2. 

We assume henceforth that H is sufficiently small that 

4rh\\A\\ 2 < 1/2. (3.5) 

The estimate {\A\y) 2 < \\A\\ 2 y 2 implies < ||v4|| 2 y 2 , for k = 1, 2, • • ■ , d. From this, 

we conclude 



(h\\A 



\2\ 



I < V 2 |j;j,'^ /2 J V 2lml exp(-y 2 /2) \Hi(A; y)\ 2 dy. (3.6) 
We next need an estimate on \Hj(A\ y)\. By Remark 5 after the definition in Section 2, 

(-A + y 2 ) KM; y)e~ y2 ' 2 = m+d) HM; y)^' 2 . (3.7) 

This equation states that 7ij(A; y) e~ y2 ^ 2 is a eigenfunction of —A + y 2 corresponding to the 
eigenvalue 2\j\ + d. Introducing normalization factors, we conclude that 

2-W 2 (j\)- 1 / 2 n- d / 4 H j (A;y)e~ yV2 = E h VkMvM ■■■ ¥>k d (Vd), (3-8) 

l*l=lil 

2 d 2 
where tpkiy) = 2 _fc//2 (A;!) -1 / 2 7r -1 / 4 Hdy) e~ y I 2 is the normalized eigenfunction of — -7— +y 2 

ay 1 

corresponding to the eigenvalue 2k + 1, and the coefficients b k satisfy J2\k\=\j\ \b k \ 2 = 1- We 
can thus deduce an estimate of Tlj(A; y) from an estimate for usual Hermite polynomials 
H k (y): 

\H 3 (A;y)\ 2 = I E hH kl ( yi )H k2 (y 2 )---H kd (y d )\ 2 

\k\=\j\ 

< J2 \H kl ( yi )H k2 (y 2 )...H kd (y d )\ 2 . (3.9) 
l*l=b"l 

The Hermite polynomials, in turn, satisfy the following bounds. 



15 



Lemma 3.1 If \y\ > \/2k + 1, then 

\H k (y)\ < 2 k \y\ k , (3.10) 
and there exists a numerical constant k > 0, such that for all y 6 M, 

\H k (y)\ < K2 k ' 2 Vk\ey 2 ' 2 . (3.11) 

Proof: The second statement is well known. See formula 22.14.17 or [|j], formula 
8.954.2. 

By symmetry, it is sufficient to prove ( |3.10D for y > 0. It is well known that the k th 



eigenfuction of the harmonic oscillator is non-zero in the classically forbidden region \y\ > 
\/2k + 1. It is also well known that H k (y) = 2 k y k + Pk{y) where Pk{y) is a polynomial of 
degree k — 2. Thus, we have 

ff*(y)/(2Y) = 1 + B k (y), 

where B k satisfies 

B k (y) = 0(1 /y 2 ) for large y, 

and 

B k {y) > -1 if y > V2k + 1. 
Using standard properties of the Hermite polynomials, we see that 

B' h {y) = (yH>(y)-kH k (y))/(2 k y k+1 ) 

= (2kyH k _ 1 (y) - k H k (y))/(2 k y k+1 ) 

= k(2yH k ^(y) - (2yH k ^(y) - 2 (k - 1) H k _ 2 (y)))/(2 k y k+1 ) 
= k(k-l)H k ^(y)/(2 k ^y k+1 ). 



So, if y > V2F+T, we see that B' k (y) > 0. Thus, y > \/2k + 1 implies -1 < B k (y) < 0, 
which implies fl3.10|) . I 



In ( p.6| ), we use the multinomial expansion 

|n| = |m| 

to obtain 



Z '"^- 7r/ |„|=| ro | |fc|=|j| V ' ll ' Ld/ r=l J 

Using Lemma \5.1[ the one dimensional integrals in the final factor here can then be estimated 
as follows (without indices): 

J y 2n \H k (y)\ 2 e W (-y 2 /2)dy 
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< f n 2 2 k k\ y 2n exp(y 2 /2) dy + ( A k y 2(n+k) exp(-y 2 12) dy 

Jy 2 <2k+l J 

< K 2 2 fc A , ! _L_( 2A , + ^n+l/2 exp ^ k + 1/2) + 4k 2 (k+n+l/2) J z 2(k+n) e -z 2 dz 

_ 2 k+1 2 fc+l/2, ,/0, , , ,fco(fc+n+l/2) 1/2 ( 2 ( k + n ) ~ 1 ) ! 

= _^_ K 2 efc+ l/ 2A;!(2A; + ^n+l/2 + 2fc -n + 3/2^1/ 2 (2(fc+^-l)! ^ 

Here we have used 



z 2( k+ n) -# dz = n X/2 



l-3-5---(2(Jfe + n) - 1) 



2 fc+n 

Now, Stirling's formula guarantees the existence of a > 0, such that 

a"n n < 7i! <n n , 
for all integers n > 1. Using this in ( p.!2|) , we obtain 

y 2n \H k (y)\ 2 exp(-y 2 /2) dy < f f k+n (k + n)\, 

for some constants /o and fx. From this we conclude (restoring the indices) 



(ft||A|| 2 ) H v v fdf H+b| / |m, \ m .... w 

J M = H |fe|=bl 

In this expression, we have 

-i 



-i 



- ^' + bl)!(, i + ni H + bl kd + nd 

< (\m\ + \j\)\ 

We further use £ ( ' m ' ) = rf |m| and ( P3[) in flgTTJ) to obtain 

|n|=[m| "' ^ 

; < (f,|WP)H ( AV (/ld) H f^Y" 



< (ft||A|| 2 )H f A) (2fxdy m \(e d fx)^ ( . |j| . ) \m\\ 



Ji ■■■ Jd 



= (h\\A\\ 2 d)M ft* & 



' 1 ji ■■■ jd 



m 
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where fa, f 3 , fi are numerical constants 
Estimates (|3.2| ) and (|3.16|) imply 

D m V(((x,a)) 



ml 



(x — a) m ipj(A, B, h, a, i], x) 



3d 



So, by the Schwarz inequality, ( |1.4j ), (|3.9|), and ( |3.18| ) 



D m V(({x,a)) 



ml 



x — a) m ipi(x, t, h) 



<EE 

n=0 \j\=n 



D m V(C(x,a)) 



ml 



[x — a) m (fj(A, B, h, a, i], x) 



< M 2 exp(8T(6 2 d + a 2 ))(h\\A\\ 2 (d/6 2 )f 3 )^f^\m\\ E E /l 

n=0 \j\=n 

7(0 

= M 2 ex V (8T(5 2 d + a 2 ))(h\\A\\ 2 (d/5 2 )f 3 )^f*\m\\ E ftd? 

n=0 

= M 2 exp(8r(6 2 d + a 2 )) (h\\A\\ 2 (d/6 2 )f 3 )^ f d 2 \m\\ f 

\ J4C1 — 1 

< M 2 exp(8T(5 2 d + a 2 ))(n\\A\\ 2 (d/5 2 )0 m \f^\m\\2(Ud)^\ 



3d 



(3.17) 



The last step depends on the following fact: We can assume without loss that q = f^d > 2, 
so that for any positive integer p, we have —1 < q p+1 — 2q p , and hence, q p+1 — 1 < 2q p+1 — 2q p . 
Thus, 

(q p+1 - l)/{q - 1) < 2q p . (3.18) 

We use the hypothesis \m\ =1 + 2 and note that |m|! < (/ + 2) l+2 . We then conclude 
that there exist constants go and g±, that depend only on d and J, such that (|2.13|) holds. 
This implies the lemma. I 



4 Proofs of Lemmas 2.3 and 2.4 

The proof of Lemma 2.3 relies on two preliminary lemmas. 



Lemma 4.1 The matrix elements of {x — a) m satisfy 

| ( ifj(A, B, h, a, i], x), (x- a) m ip k (A, B, h, a, r), x))\ 

< ftM/ 2 (V2d)H \\A\\™ J(\h\ + l)(\k\ +2).-. (\k\ + |m|), 



(4.1) 



and 



( <fj(A, B, h, a, 77, x), (a? - a) m (f k (A, B, h, a, 77, x) ) = 0, if 



\J\ ~ \k\ 



> m 



Proof: For i — 1, 2, • • • , d, we can use equation (3.28) of J7| to express (xj — aj) in terms of 
raising and lowering operators. Doing so, we obtain 



(xi - Oi) ip h (A, B, h, a, 77, x) 



Ap\/k p + l(pk'(p)(A B, h, a, 77, x) 
P =i 

d - I— 

+ yh/2Y^ A ip ykp(pk»(p)(A, B, h, a, rj, x), 

P =i 



where k'(p) and k"(p) are constructed from k by replacing the component k p by k p + 1 
and k p — 1, respectively. Thus, (x^ — ai){p k (A, B, h, a, 77, x) may be written as a sum 
of 2d terms, each of which has the form b q ip q (A, B, h, a, 77, x), where |g| < \k\ + 1 and 



\b g \ < y/h/2\\A\\ y/\k\ + 1. 

From this and a simple induction, (x — a) m (fk{A, B, H, a, 77, x) may be written as a sum 
of (2(f)!" 1 ! terms, each of which has the form b q (p q (A, B, h, a, 77, x), where \k\ — \m\ < \q\ < 
\k\ + \m\ and \b q \ < (h/2)\ m \ \\A\\\ m \ yj(\k\ + l)(|fc| + 2) • • • (jjfej + |m|). 

This implies the lemma. I 



In the next lemma we use the shorthand X m (t) to denote X m (l, t, h) of Lemma pO 



Lemma 4.2 For each N > 0, there exists D(N) < oo, such that 
sup \{X m ^{t l )X m >{t 2 ) ■■■ X m "[t q )) jk \ 

{j,k: |j|,[*|<J(0} 

C\ |mi|+|ma|— +|m 9 | 
D(N)Vh SUp \\A(t)\\) J(J)(\mi\+\rm\+-+\m q \)/2 i 
te{ti,t a ,-t q } J 

for any collection mi, 777.2, ■ ■ ■ , m q of multi-indices that satisfy \rrij\/J(l) < N. 



(4.2) 



Proof: The matrix X m is obtained from X m by restricting its indices to \k\ < J (I). 



Using Lemma 4.1, we see that for such values of the indices, 



xf k {t) 



< 



^d P(t)||)H j(/)H/2 v /(i + l/J(/))(l + 2/J(/))).-.(l + |m|/J(0). 
Thus, \m\/J(l) < N implies 

Xf k {t) < (V2hd\\A(t)\\) lml (I + N) lml/2 J(l) lml/2 . 
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We estimate the absolute value of the j, k matrix element of the product 
\{X^{t x )X m \t 2 ) ■■■X m «(t q )) jk \ 



< 



E 



Xf^tx) x™ 2 P2 (t 2 ) ■ • ■ x™ 9 _ lk (t q ) 



\Pi\ < J{1), t = l,---,g-l 
\pi -j\ < \mi\,\p2 -pi\ < Imal,-"! \k-p q -i\ < \m q \ 



(4.3) 



by bounding the absolute values of each of the matrix elements X™ l _ ip .(ti) by 

X™l in (ti)\ < {\Pmd\\A{U)\\)\ mi \ (l + iV)l m 'l/ 2 J(/)l m 'l/ 2 , (4.4) 

and multiplying by a bound on the number of terms. 

To estimate the number of terms, we first note that for any multi-index r, the number of 



multi-indices indices p that satisfy 



\p\ ~ \ r \ 



< \m\ is equal to the number of vectors v with 
m\. This is bounded by the number of vectors 



integer components, such that | X^=i v i\ < 
with integer components, such that \vi\ < \m\ for % = 1, 2, • • • , d, which is (2\m\ + l) d . 
Thus, the number of terms is bounded by 



E 



i 



\pi\ < JQ), i = l,---,9-l 
\P1 ~j\ < |mi[,|p2 —Pl\ < \m 2 \,--- , \k — Pq-l\ <\in q \ 

< {2\m x \ + l) d (2|m 2 | + l) d ■ ■ ■ (2\m q \ + l) d 

< ^ e 2-jd(|mi| + |m 2 |H h|m g |) 

The final inequality follows because 1 + 2\rrij\ < e 2 ' mj '. 

We obtain the lemma by using (|4.4| ) and (|4.5|) to bound ( |4.3|) . 



(4.5) 



Proof of Lemma |2.3| We mimic the proof of Lemma [4.2| . Let m = \mo\ + \mi\ + |m 2 | + 
. . . + \m q \. By Lemmas [4.1| and 42 , 



[Vh sup \\A(t)\\ 
V te{t ,ti,-,t q } 



(X^i^X^^X^ih) ■•• X m «(t q )) jk 



Vh sup \\A(t)\\) ]T \x™°(t )\ \(X m -(t 1 )X m *(t 2 ) ■■■ x m «(t q )) rk 

te{to,tl,™.*g} / | r |<j a) 



< E 2l-l/ 2 ^1 + 1)^1 + 2) ■■■(|j| + |m |) 

x £)(jV)' mi ' + ' m2 ' + '" + ' m<J J(/)(l mi l + l m2 l + ---+l m 9l)/ 2 



r < J{1) 
\j -r\< \m \ 



< E 2 ' mol/2 Al) lmol/2 V(! + + 2/J(0) •••(! + \m \/J(l)) 



\r\J(l) 
\j -r\< |m | 



20 



x £)(jV)' mi ' + ' m2 ' + '" + ' m9 J(/)(l mi l + l m2 l + -" + l m 9l)/ 2 



< ^ 2 |mo|/2 (1 + jV) |m,)l/2 D(N) 1 



m 1 \ + \rri2\ + ...+\m q \ JH\ (\ m l [+|"ia|+—+|m g |)/2 



|r| < J(l) 
\j -r\< |m | 



< 2l m ol/2 ^ _|_ jyymol/2 £^^jy^|m l |+|ma|+...+|m,| j^y\m 1 \ + \m 2 \+...+\m q \)/2 

b'-r]<|m | 

Since the sum in the last line is (l + 2|mo|) d < e M ' m °', we obtain the desired estimate. 



We now turn to the proof of Lemma |2.4| . Our first step is to study the combinatorial 
factor F p (n, q) in the one dimensional case. 



Lemma 4.3 Let rrij denote positive integers, and let 

G p (n,q) = £ 1. (4.6) 

1 < mi , m,2 , ■ ■ ■ , m q < p 
mi + m% H + m 9 = n 

27ms quantity is zero if n < q or n > qp. Otherwise, it satisfies 

f ( n ~]) if q<n<[q(p + l)/2] 

(( 9lP g _! ) ^ [q( P +l)/2} + l<n<qp 

Proof: If n < q or n > qp, the sum in ( |4.6| ) contains no terms, so G p (n, q) = 0. 

The quantity G p (n, q) is the number of ways that the number n can be decomposed as n = 
mi + ni2 + . . . + m q with each rrij satisfying 1 < nij < p. We uniquely associate to each such 
decomposition of n, a corresponding decomposition of n' = q{p+ l)—n = m\ +m' 2 + . . . +m' q 
by setting m'- = p + 1 — rrij. This association is a one-to-one correspondence, so we see that 
G p (n, q) satisfies the symmetry relation 

G p (n,q) = G p (q(p + l)-n,q). (4.8) 

As a consequence, the second inequality in ( fL7|) (with [q(p + l)/2] + l<n< qp) follows 
from the first (with q < n < [q(p + l)/2]). 

To prove the first inequality in (|4.7|) we drop the condition rrij < p and exactly calculate 
the resulting function. Dropping the upper bound on the rrij, we clearly have G p (n, q) < 
G(n, q), where 

G(n, q)= J2 L 

1 < mi, m2 , m q 
mi + rri2 + ■ ■ • + m q = n 

So, the lemma will be proved once we establish 

G(n, q) = "J) . (4.9) 
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To prove this, we use induction. Formula ( |4.9| ) is trivial to verify for all n > 1 when 
q — 1 and for all n > 2 when q = 2. Assume (}4.9| ) has been verified whenever n > g — 1, and 
suppose n> q. We have 



n—q+l n—q+l 

G(n, q) = G(n-mi,q-l)= E 

mi=l mi=l 

so the induction step will be complete once we show 

n—q+l 

E 

mi=l 



n — rrii — 1 
q-2 



( n — mi — 1 \ 




/ n - 




I g-2 J 




\q- 


-I) 



(4.10) 



To prove ( |4.10|) , we do an another induction on n > q. For n = q we trivially have 



Now assume ( [4.10D is true for some n> q. For the case of n + 1, we have 



(n+l)- g +l 

E 

m=l 



(n + 1) — m — 1 
g-2 



n—q+l 

E 

m=l 
n— g+1 

E 

m=l 
n— g+1 

E 

m=l 



n — m 
q-2 

n — m 
q-2 

n — m 
q-2 



+ 



+ 



' n — m — 1 
9-2 



E 

m'=0 



> + l)-{(n + l)-g+l}-l 
q-2 



q-2 
q-2 



+ 1 



n — 1 
g-2 

n — 1 
q-2 

n 
q-1 



E 

m'=l 

n — 1 
q-1 



( n — ml — 1 \ 




( 9 




I 9-2 J 









+ 1 



This proves ( 4.10|) and completes the proof of the lemma. I 

Our next step is to obtain an estimate for G p (n, q) that depends only on n. 
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Lemma 4.4 Let G p (n, q) be defined by (^..dj) . There exists a constant C 7 , such that p > 1, 
q > 1, and n > 1 imply 

G p (n,q))<(C 7 r. (4.11) 

Proof: We note that G p (n, q) is zero unless q < n < qp. Furthermore, G p (q, q) = 
G p (qp,q) = G2(2,l) = 1 and Gi(n,q) = 5 nq . So, we need only prove existence of C7, 
such that p > 2, q > 1, (p, q) ^ (2, 1) imply 

sup (G^q)) 1 ^ < C 7 . (4.12) 

q<n<pq 

To prove this, we first study the case where q < n < [q(p + l)/2]. Since q < n, we have 

q l l n = exp(ln(g)/n) < exp(ln(n)/n) < e 1 ' 6 . 
Thus, by Q4.7|) , the condition q < n < [q{p + l)/2], and ( |3.13|) , we have 



< q ^ ( ("-D' Y" 

\ql(n- q)\ J 

< eVe / / 71 m TT- ( 4 - 13 ) 

An explicit computation shows that 

d 71 ft q 

On qil n in - g) ( n ~i)/ n ~ qil n (n - q) n 2 ^ 9 ' 7 ' 

From this, we deduce that the right hand side of (|4 . 1 3|) attains its maximum at n — 2q. 
Evaluating that maximum, we obtain 

sup (G p (n, q)) 1/n < -e 1/e . (4.14) 

q<n<[q(p+X)/2] a 

For [q(p + l)/2] < n < pq, the number n' = q(p + 1) — n satisifes q < n' < [q(p + l)/2] 
and n'/n < 1. Thus, by ( |4.8| ) and ( [4. 14] ), we obtain 



(G p (n, q)f' n = {G p {n\q)fl n 
< I - e 1/e 



a 

Thus, 

sup (G p (n,q)) 1/n < -e 1/e . (4.15) 

[<j(p+l)/2]<n<p<j a 
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Inequalities (|4.14 ) and ( |4.15| ) imply the existence of C-j for which ( |4.12| ) holds. This 
implies (|4.11|) ,and the lemma is proved. I 

We now generalize Lemma [4.4| to the multi-dimensional case. 



Lemma 4.5 Let F p (n } q) be defined by ( \2.1o\ ). For all p > 1, q > 1, and n > 1, we have 

F p (n,q) < C^ d \ (4.16) 
where C-j is the constant of Lemma \^-4- 



Proof: With rrij temporarily denoting numbers instead of multi-indices, we define 



TJn, q) 



E 



< mi, 7712 , • • ■ , iriq < p 
mi + ni2 H h m q = n 



which is zero unless < n < pq. 

By defining m'- = nij + 1, we see that every decomposition of n as n = m 1 + m2 + • • ■ + m q , 
with < nij < p, corresponds uniquely to a decomposition ofn + qa.sn + q = m' 1 + m' 2 + 
• • • + m' q , with 1 < m'j < p + 1. Therefore, we have the identity 



r p (n, q) = G p+1 (n + q,q). 



(4.17) 



We now let the rrij denote multi- indices and let rrij(k) denote the k th component of rrij. 
Then, using (|4.17|) and Lemma 4.4 , we easily obtain 



F p {n, q) < 1 = r p(™> <l d ) ^ C 

< rrij (k) < p 

ELi ELi m ^ k ) = n 



(qd+n) 
7 



Proof of Lemma |2.4| : Suppose 2 < q < I + 2, and let L q (l) satisfy ( |2.17| ), i.e., 



LS) < E (C 2 hir/ 2 F l+1 (n,q). 

n Q- n=n-2 



By Lemma [4.5| and (|3.13| ), there exist Cg and Cg, such that 
(dC 7 t 9(m) 



^(0 < 



y {c 2 ni) n/2 c? 



< 



g g(Z+l) 



(M 9 n = 



E (c 8 ^ 1/2 / 1/2 ) n . 



:3«-2 



(4.1* 
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Note that we can take 



C 8 = rnax{l, yJC 2 C 7 }, (4.19) 
C 9 = max{l, C X C$}, (4.20) 

because we can assume C§ > 1 and Cg > 1 without loss of generality. 

We arbitrarily choose two positive numbers a < 1 and ft <1. We then define 

5 = ci min i a ' cfcf^r (4 - 21) 

We henceforth assume I = 1(h) is chosen to satisfy (|2.18 ) with this value of g. With this 
choice, we have 

C s h 1/2 l 1/2 < a. (4.22) 
By summing a geometric series in (|4.18|) , we see that 

- (1 -g£ci)» ( C8ftl/2 ' 1/2 ) 3 " 2 - (4 - 23) 

As a function of g, 

(g) = exp ( gln (V9)) 
is maximized at g = Z/e, where it has the value e l ^ e . Thus, 

©' * " 

This, g < Z + 2, and ZiZCf < a 2 < 1 imply 

1 1 flV e l ' e 

< 



From this, Cg > 1, and Cg > 1, we conclude 

r l+2 r 2(l+2) i/ e 

L q (l) < ° 9 ,f 8 . 6 (C 8 h^l^) 1 - 6 
[I — a) K 7 

- v e c ^ n 1 ) (i -a) (hirer 

By ( 14.211) , we have 

eV'CgClhWl 1 ' 2 < (3 < 1. 



9 1 < Z < I, 



Since 

h 

25 



(4.24) 



we see from ( 4.24|) that if h < g/2, then we have 
L q (l) < e-l^l' 



r<2 

°9 



< 



2 3 eC 2 



(l-a)g*CI 



(1-a) (g-hfCi 

-\\n{p)\g/h 



(4.25) 



Note that the factor 2 3 can be replaced by (1 + e), with e arbitrarily small, by taking h 
sufficiently small. M 

3 1 B 2 

Remark We can choose a and 8 to satisfy —-r < a < —rr- Then, since — = — < 

eV e - e^h ClCle 2 / e 

e 2/e 



< a , we obtain the conclusion of Lemma |2^J with 



9 
7i 
C 3 



/? 2 



C 9 2 C 8 6 e 2 / e ' 

| ln(/?)| /? 2 

C 9 2 C 8 6 e 2 / e 

2 3 e 1+6 / e C 9 8 C 8 10 
(l-e-Ve)/36 ■■ 



(4.26) 
(4.27) 
(4.28) 



where C$ and Cg are related to C\, C2, and C7 by ( |4.19| ) and ( |4.20| ). Note that C-j is purely 
combinatorial and has no time dependence. 



5 Proof of Theorem 1.2 



To prove Theorem |1.2|, we revisit the proof of Theorem [LT] and make the T dependence 
of all constants explicit. We then allow T to grow with h with the restriction that our 
approximation remain close to the actual solution as h —>■ 0. 



Lemmas 4.3, 4.4, and 4.5 have no time dependence, and the time dependence of Lemmas 



2A\ , |2.2| , |2.3| , |3.1| , |4.1| , and |4.2| has been made explicit in their conclusions. 

Since V is bounded below and energy is conserved, \a(t)\ grows at most linearly with 
time. Thus, there exist v i > and v > 0, such that v(t) defined by ( p.33|) satisifies 



v(t) < Vq exp(t>xt). 
From assumption ( |1.13| ), it follows that 

w(T) = sup \\A(t)\\ < N exp(AT). 

0<\t\<T 



(5.1) 



(5.2) 



In the sequel, we denote all innessential constants that do not depend on T by the same 
symbol c. 
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Consider Lemma |2T2|. For z G Cg(C(x,a(t))), we can prove the existence of p > such 



that 



D m V{C{x,a{t))) 



ml 



[x-a{t)) m ^{x, t, h) 



< c exp{p\t\) \cy/h(l + 2) exp(A|t|) 

Indeed, with our bound on [ V(^) | we have 

D m V(((x,a(t))) 



1+2 



(5.3) 



ml 



< c exp(p\t\) exp(c|x — a(t)\) 



instead of (|3.2|) . With this and the estimate exp(cy) < exp(q/ 2 ) + exp(c) for all c > and 
y > 0, we see that the estimates in the proof of Lemma |2]2| are still valid. This yields (|5.3| ). 
Consequently, the constants C4 and C5 appearing in (|2.24j) satisfy 

C 4 (T) = c exp(pT) 
C 5 (T) = cexp(AT), 

where p > 0. By using estimates (|5.1|) and ( |5.2j) in ( |2.34j ) and ( p.35| ), we see that when we 
apply Lemma |2]4|, the constants C\ and C 2 satisfy 

d(T) = cexp(pT) 
C 2 (T) = cexp(AT). 

We still must determine the T dependence of C and 7 of Theorem |1 . 1| . We do this by 
determining the T dependence of g,7i and C3 in Lemma [2.4| . Equations ( 4.19 ) and ( 4.20 ) 
yield C§ and C9, in terms of which the above constants are determined. We find that 

C 8 (T) = cexp(AT/2) 

C 9 (T) = cexp(pT). 

Consequently, from flPSD , QOTD , (|428|) , we have 

(7(T) = c exp(— i/T), f = 2p + 3A 

7i(T) = cexp(-z/T), 

Cs(T) = c exp(/iT), p = 8p + 5 A. 

To arrive at these conclusions, we imposed various conditions. Those conditions were 
( P-5|) , that now requires h < cexp(— 2 XT), 
( j2.36| ), than now requires cexp(2AT) exp(— vT) < 1, and 
h < g(T)/2 used in (f4.25|) , that now requires h < cexp(— vT). 
These are all satisfied, provided we take 

T(h) < T'|ln(^)|, 

with T' > sufficiently small. Using this in our estimate for the error term £(t, h), we see 
that for sufficiently small T' , \t\ < T' \ hx(h)\ implies 

\\^{x,t,h)-*(x,t,n)\\ L . m < c"exp(-y/r), 

for some 7' > 0, o > 0, and C > 0. I 
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